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A single impurity problem is investigated for multiband s-wave superconductors with 
different sign order parameters (s±-wave superconductors) suggested in Fe-pnictide super- 
conductors. Not only intraband but also interband scattering is considered at the impurity. 
The latter gives rise to impurity-induced local boundstates close to the impurity. We present 
an exact form of the energy of the local boundstates as a function of strength of the two 
types of impurity scattering. The essential role of the impurity is unchanged in finite number 
of impurities. The main conclusions for a single impurity problem help us understand effects 
of dense impurities in the s±-wave superconductors. Local density of states around the sin- 
gle impurity is also investigated. We suggest impurity site nuclear magnetic resonance as a 
suitable experiment to probe the local boundstates that is peculiar to the s±-wave state. We 
find that the s±-wave model is mapped to a chiral d x i_ y i ± id xy -w&vc, reflecting the uncon- 
ventional nature of the sign reversing order parameter. For a quantum magnetic impurity, 
interband scattering destabilizes the Kondo singlet. 

KEYWORDS: multiband superconductivity, impurity, boundstate, Fe-pnictide superconductors, 
s±-wave superconductivity, nuclear magnetic resonance, numerical renormaliza- 
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1. Introduction 

The investigation of impurity effects on superconductivity has been developing for a long 
time. In case of conventional s-wave (BCS) superconductors, Anderson showed that nonmag- 
netic impurities change neither the superconducting transition temperature (T c ) nor the gap 
of the superconductor. 1 ) It is first pointed out by Abrikosov and Gor'kov that the magnetic 
impurities cause gapless behavior. 2 ) They reduce the superconducting energy gap and suppress 
T c . As a result, they also give rise to finite density of states inside the superconducting energy 
gap. 3-5 ) In the same manner as the dense magnetic impurities, a single impurity brings about 
localized boundstates inside the energy gap. 6 ' 7 ) While the problem of a classical spin can be 
solved exactly, 5 ' 8 ' 9 - 1 a quantum spin involves us in a many-body problem of the Kondo effect. 10 ) 
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The latter case had been studied by various theoretical methods 11-13 ) and was finally solved 
for the conventional s-wave superconductivity using the Wilson's numerical renormalization 
group (NRG) method. 14 ^ 16 ) 

In contrast to the BCS superconductors, nonmagnetic impurities destroy unconventional 
superconductivity. For instance, Zn impurities in cupper oxide high temperature supercon- 
ductors induce additional finite density of states inside the superconducting energy gap, which 
accounts for the temperature dependence of nuclear magnetic relaxation (NMR) rate. 17 ) A 
single impurity problem in d x 2_ y 2-wave superconductors was also studied. It was found that 
low energy states appear with four- fold symmetry near the impurity. 18 ~ 22 ) About a single mag- 
netic impurity, we previously investigated a quantum spin in unconventional superconductors 
using the NRG method and focused on a fully gapped chiral superconductor expressed by 
p x ± ipy-wave or d x 2_ y 2 ± id^-wave type order parameters, 23 ~ 27 ) where orbital effect of the 
Cooper pair plays an important role. 

Recently, Kamihara and coworkers discovered a new Fe-pnictide superconductors. 28,29 ) It 
is suggested theoretically that antiferromagnetic spin-fluctuations arising from the interband 
nesting favor a multiband s-wave superconductivity with different sign order parameters 30,31 ) 
that is called s-t-wave here. For this new type of multiband superconductivity, interband scat- 
tering is important. It affects NMR relaxation rate, 32,33 ) can suppress T c , 34,35 ) and generates 
impurity-induced states inside the energy gap 36-38 ) similarly to a magnetic impurity in BCS 
superconductors. 

The purpose of this paper is to investigate the single impurity problem to understand the 
novel properties of the s±-wave state. For the single impurity, we can obtain an exact solution 
that helps us understand the properties of many impurity case. The following points will be 
clarified in this paper: (1) An explicit form of energy of the impurity-induced boundstates 
is presented as a function of strength of the interband and intraband scatterings. Spatial 
dependence of the local density of states is shown around the impurity. (2) The pair breaking 
effect of the interband scattering is interpreted from an effective single band model. Relation to 
chiral d x 2_ y 2 ± id^-wave superconductors is also discussed. (3) Quantum effect of a magnetic 
impurity in the s±-wave superconductors is analyzed by the NRG method. 

This paper is organized as follows. In §2, we study nonmagnetic impurity and discuss 
effects of interband and intraband scatterings on appearance of the localized boundstates. In 
§3, we focus on an identical multiband case and discuss the same problem from a point of 
view of an effective single band model. Then our theory is extended to a quantum magnetic 
impurity in §4. The last section gives summary of our results. We assume h = 1 and &b = 1 
throughout this paper. 
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2. Nonmagnetic Impurity 

2.1 Formulation 

Let us begin with the following Hamiltonian of a continuum model for multiband super- 
conductivity with impurity scatterings: 



H » = E/ dr^ CT (r)e M (-iV)^ CT (r) - A M J dr [v-J T (r)^(r) +^(r)V„ T (r) 
H ' = E E/ dr^(r)^ )<r<r ,(r)^(r). (1) 

fXfx'=± era' 

Here, W M is the BCS Hamiltonians for the //(= ±) band, in which tpfxcr{ r ) is the field operator 
of the conduction election of cr(=j, j) spin and e M (iV) = — V 2 /2m M — E'p is the operator of 
the kinetic energy for the \i band measured from the Fermi energy where represents the 
band dependent mass of the conduction electron. A M is the fi band superconducting order 
parameter. We assume that A M is a real value and that the sign of the order parameters are 
different (A + A_ < 0) between the two bands. H' represents the Hamiltonian for the impurity 
scatterings. U^i )Crrj i(r) is the amplitude of the scattering between the fi band electron with a 
spin and the fj,' band electron with a' spin. A single nonmagnetic impurity is located at the 
origin of the coordinate. The scattering amplitude is given by 

U^^(r) = U^5 aa ,8(r). (2) 

Here, 5 aa ' and S(r) are the Kronecker delta and Dirac delta functions, respectively. The 
fi = \j! components are for the intraband scattering, while the fi ^ fi' components are for 

the interband scattering. We assume that U^i is a real value and U-\ = U |_. For the 

nonmagnetic impurity, we define the following thermal Green's function in a 4 x 4 matrix 
form: 

G(T,r,r') = -(T T *(r,T)&(r')), (3) 
where \&(r) and *ff\r) are 4 dimensional vectors. The latter is defined as 

*V) = ( V^M <Mr) Vl T (0 ^(r) ) . (4) 

The imaginary-time Heisenberg representation is defined by 

*(r,r) = e nT ^{r)e- HT . (5) 

In the absence of the impurity scattering, the unperturbed Green's function in the Fourier 
transformed form is given by 

n r n ( (: *^'- k ) \ m 
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Here, uji = 2nT(l + 1/2) is the Matsubara frequency for fermion. G± is unperturbed 2x2 
Green's functions for the ± band. It is given by 

iwi + e± )k p 3 - A±pi 



G±(m,k) 



oof + e 2 , . + Al 



(7) 



where p a (a = 1, 2, 3) are the Pauli matrices for the particle-hole space. e±^ = k 2 /2m± — Ep 
is the band dependent kinetic energy. The real space Green's function is obtained as 

G (m,r,r') = i^e ifc (— ')<?„(*;,, fc) 



(8) 



Here, represents the system volume. 

In the presence of the impurity scattering, the Green's function is calculated exactly as 19 ) 



-l 



G(kj h r, r') = G (i^, r, r') + G (iwj , r, 0)17 |i - G (m, 0, 0)C/j G (iw,, 0, r'). (9) 
Here, £/ is given by 



17 



V 








U+- 








-U++ 





-u+ 







L7__ 








-*7+- 





—[7- 



(10) 



Go(iw/,0, 0) in eq. (9) is calculated as 

G o (iw,,0,0) = ^ J]Go(iwi,fc) 



G±(iwj,0,0) =-7TiV ± 



G+(iw,,0,0) 

G_(iw,,0,0) 



iw; - A ± pi 



(11) 



V^ 2 + A ± 

Here, N± represents the density of states per volume at the Fermi energy for the ± band, 
respectively. Go(icu;,r,0) and Go(iu^,0,r) in eq. (9) are calculated as 

/ G+(iw,,r,0) \ 



G o (m,r,0) = G (iuJi,0,r 



where 



V G_(iw,,r,0) 

G±(iiu l ,r,0) = I^e ifc -G±(i^,fc). 



(12) 



(13) 



For isotropic two dimensional conduction electron systems, it is expressed by Bessel functions 
as in eq. (A-6). Details of the integration are given in the Appendix. 

The local density of states at position r is given by the Green's function. Since we consider 
nonmagnetic scatterings here, density of states are same for a =|, [ spins. They are expressed 
as 



Ni(E, r) = N[(E, r) = N intia {E, r) + N intcI (E, r), 



(14) 
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N intia (E,r) = a u (E,r) + a 33 (E,r), N intcI {E,r) = a 13 (E,r) + a 31 (E,r). 
Here, -/Vmtra (-/Vinter) is the intraband (interband) contribution. a m n(E,r) is defined by 

a mn {E,r) = Im G{\uji ->■ E + id, r, r) , (15) 

7T L J mn 

where 5 is a positive infinitesimal small number and [• • • ] mn represents the mn matrix element. 
For dense impurities, the local density of states is uniform after averaging over the impurity 
positions. In this case, only the intraband contribution remains, since the a± 3 and 031 terms 
vanish after integrating over the coordinate due to the orthogonality of the wavefunctions 
of different conduction bands. For the single impurity, however, these terms remain in the 
presence of the interband scattering. 

2.2 Impurity-induced local boundstates 

The interband scattering connects the two bands with different sign superconducting order 
parameters. We can expect boundstates as in case of many impurities. 36 ) Since an important 
result is not altered by details of the multiband structures, we study here a case of |A + | = 
|A_| = A. In this case, we can express energy of the boundstates explicitly. This gives us 
useful information about the density of states at low energy when many impurities are taken 
into account. 

Energy of the impurity-induced local boundstates is determined by poles of the Green's 
function given in eq. (9). Solving |1 — Gq{\uji — ► E, 0, 0)U\ = 0, we can determine the bound- 
state energy positions. They are expressed explicitly as 



E - t A. 1 + ^ + "-- + ^- n++ "")'- 2 "+-. (16) 
\| 1 + u\ + + u 2 _ _ + {u\_ - n ++ M__) + 2u\_ 

Here, the sign of the energy (±) corresponds to particle-like and hole-like excitations, respec- 
tively. n ++ , u , and n_| are defined by 

u++ = 7TN+U++, u__ = vriV_L7__, u+_ = ir y / N + N-U + -. (17) 

Let us discuss various cases below. 

2.2.1 Effect of interband scattering 

First, we restrict ourselves to effect of the interband scattering and discuss an identical 

multiband case here (N + = iV_, u ++ = u , and £ + = £_). When we put u ++ = u = in 

eq. (16), we obtain 

1 — u 2 

E = ±Eb, E b = -A +- sgn(^ + _). (18) 

There are two boundstates corresponding to particle-like and hole-like excitations. They in- 
tersect at u+_ = 1 (or at E = 0) as shown in Fig. 1(a). In the unitary limit (u + _ — > 00), 
the boundstate energies come close to the superconducting energy gap. This means that there 
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Fig. 1. (Color online) (a) Interband scattering dependence of the boundstate energy for u ++ = 
= 0. Solid and dashed lines are for E = ±£7bj respectively, (b) Local density of states at 

the impurity site (r = 0) described by eq. (19) for u_| = 0.5. We introduced a finite broadening 

factor (r = 0.001 A). 



is no boundstate in the unitary limit. We notice that eq. (18) has the same form for a sin- 
gle classical spin in conventional s-wave superconductors. 5,9 ) In this sense, the nonmagnetic 
interband scattering in s-t-wave superconductors plays the same role of a classical spin in 
conventional s-wave superconductors. We will discuss this point in §3. 

Let us see the local density of states at the impurity site (r = 0). For the identical 
multiband, it is expressed as 

N a (E,0) _ 1 \E\VW^ 2ttA 

iV+T^T " l + u%_ E*-E% ° m ~ A) + (l + u%y {E " Eb) > (19) 
where 9{x) is the Heviside step function. The first term in eq. (19) is for continuum states, 
while the second term is for the local boundstate. Only one boundstate is visible as in Fig. 
1(b) in the identical multiband case. Another characteristic point is that the intensity of the 

density of states decreases with increase of |it_| |. This means that the wavefunction vanishes 

at the impurity for large values of \u^ |. 

We next show spatial dependence of the local density of states dividing it into the intraband 
[Fig. 2(a)] and interband [Fig. 2(b)] contributions. We can see that there are two boundstates 
inside the energy gap with Friedel oscillations for both continuum and boundstates. The 
intensity of the boundstates decay in the superconducting coherence length £. Adding both 
contributions, the density of states for the boundstate at E = E-q only remains as in Fig. 2(c) 
similarly to a single classical spin in conventional s-wave superconductors. 5 ' 9 ^ This is due to 
the character of the identical multiband. 

When the two bands are not identical, the cancelation between the two contributions is 
not perfect, which results in the two boundstates inside the energy gap as in Fig. 2(d). At the 
impurity site (r = 0), the finite value for the E = — Eb boundstate is due to the N + ^ N- 
character here. Generally, |A+| 7^ |A_|, N+ 7^ and 7^ u characters give rise to the 
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N„(E,r) N (E,r) 




Fig. 2. (a)-(c) Spatial dependence of the local density of states for the identical multiband (|A + | = 
|A_|, N + = N-, u ++ = u—, and £ + = = £) renormalized by the value in the normal state. 
The radius r is scaled by the coherence length £ defined in cq. (A-8). The ratio Ep/A = 100 [see eq. 

(A-9)] and the broadening factor T = 0.1A are used, (a) Ai ntra for u_| = 0.5 (m++ = u__ = 0). 

(b) Ajnter- (c) N a = Ai ntra + Ai n t er . (d) N a for a non-identical multiband. Set of parameters are 
chosen as |A+| = |A_|, 7V_ = 0.27V+, u ++ = u = 0, w + _ = 0.5, and £ + = 0.2^_. 



finite intensities for the two boundstates at r = 0. When those values are same for the two 
bands, there is no intensity for the E = —E-q boundstate at r = 0. The intensity at r / 0, 
however, can be finite when since the cancelation becomes imperfect there. In reality, 

the two bands are not identical and the appearance of particle-like and hole-like boundstates 
are expected inside the energy gap. It is reported that there are two boundstates inside the 
energy gap in a tight-binding model calculation. 37 ' 38 ) 

2.2.2 Effect of intraband scattering 

Next, we examine effect of intraband scattering. We consider here the identical multiband 
case. The boundstate energies are given by 

fl = +A l + u\ + -uj_ (2Q) 

^/l + (u2 + - U 2_ )2 + 2(u 2 ++u 2_ ) 

We introduce a parameter a = u ++ /u^ as the ratio of strength of the intraband and in- 

terband scatterings. Figure 3(a) shows the dependence of the boundstate energies for 
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Fig. 3. (Color online) (a) u_| dependence of the boundstate energies for various a(= u ++ /u^ ) 

values, (b) Spatial dependence of the local density of states for the identical multiband in the 
unitary limit. The ratio Ep/A — 100 and the broadening factor T = 0.1 A are used. 

various a. For a < 1, the two boundstates intersect, while they do not for a > 1. For both 
cases, there is no boundstate in the unitary limit. 

The result is different when a = 1 as discussed by Senga and Kontani. 35 ) For a = 1 
(u ++ = = u), eq. (20) becomes simple as 

E = ±E B , E B = -A—=L==sgn(u). (21) 

a/1 + (2u) 2 

In contrast to the a / 1 case, the boundstate energy becomes E B — ► for u — > oo. This 
indicates that there is a mid-gap boundstate for a = 1 in the unitary limit. We note that the 
expression of eq. (21) is same as that for a single nonmagnetic impurity in chiral supercon- 
ductors such as a p x ± ip^-wave or d x 2_ y 2 ± id^-wave type. 39 ' 40 ) For a = 1, the impurity site 
local density of states is expressed as 

N A E, 0)= I JEIVW^ ,|2,,| 
W+ + JV- l + (2»)2 E'<-El [l + (2„)2]§ 

As in the a = case, there is only one boundstate {E = E B ) in the local density of states 
due to the identical multiband character. In Fig. 3(b), we also show the spatial dependence 
of the local density of states in the unitary limit. There is a mid-gap boundstate with Friedel 
oscillations. The local density of states is much suppressed at the impurity site by the strong 
scattering. We note that both boundstates (±E B ) have finite intensities for non-identical 
multiband. 

2.2.3 Effect of non-identical multiband 

When the two bands are not identical, the following quantities are different: N + ^ 
u ++ ^ u , |A+| ^ 7^ In this case, both the particle- like and hole- like bound- 

states have a finite intensity in the local density of states as we discussed in the previous 
subsections. Besides this, energies of the boundstates change from the identical multiband 



8/23 



J. Phys. Soc. Jpn. 



Full Paper 



case as expressed in eq. (16). In |A+| = |A_| case, the boundstate energy is expressed explic- 
itly by eq. (16) as a function of u , and In both the denominator and numerator, 

the {u\_ — u ++ u ) 2 part is dominant in the unitary limit. Since the boundstate energy be- 
comes ±A, there is no boundstate inside the energy gap. In contrast to this, the boundstates 

stay inside the energy gap even in the unitary limit when u^___ = u ++ u (or U\_ = ) 

is satisfied. In this case, the boundstate energy is expressed as 



We plot the u+_ dependence of the boundstate energy in Fig. 4(a). There are two boundstates 



In |A + | 7^ |A_| case, it is difficult to express the boundstate energy explicitly. We find poles 
of the Green's function and determine the boundstate energies [see Fig. 4(b)]. Compared to the 
|A+| = | A_ | case, the boundstate energy shifts inside the smaller energy gap. We also show 
the spatial dependence of the local density of states in Fig. 4(c). There are two boundstates 
inside the smaller energy gap. Thus, the boundstate can exist even in the unitary limit when 
Uf_ = U ++ U is satisfied. 

We examined effect of a single nonmagnetic impurity with interband scattering and found 
that local boundstates appear near the impurity. When there are many impurities, the local 
boundstates overlap each other and form an impurity band as in the conventional s-wave 
superconductors with magnetic impurities. 5 ) The center of the impurity band is determined 
by the energy of the local boundstates, while the width of the impurity band is controlled by 
impurity concentration. Thus, eq. (16) is useful for knowing the in-gap state that appears in 
the density of states for dense impurities as examined by Senga and Kontani. 36 ) 

2.3 Impurity site nuclear magnetic resonance 

For the s±-wave superconductivity, boundstates appear when there is an interband scat- 
tering. Since the local boundstates exist near the impurity, we can expect those low-energy 
excitations to be detected by some local probes. Scanning tunneling microscope is one of the 
candidates. 41 ^ Besides this, we discuss here impurity site NMR 40 ) as another candidate of a 
local probe to examine the exotic superconductivity of the s±-wave state. 

The impurity site NMR relaxation rate is proportional to the following dynamical spin 
correlation function: 42 ) 



Here, S±(t, 0) is the spin operator in the Heisenberg representation at the impurity site (r = 0). 

is a retarded two body Green's function. The thermal Green's 




(23) 



in the unitary limit for u 2 



■+- = u ++ u— 




(24) 
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Fig. 4. (Color online) (a) u_| dependence of the boundstate energies for |A + | = |A_|. The solid 

(dashed) line is for u ++ = 2u__ = y/2u-\ («++ = 2u__ = ). (b) For |A_| = 0.5|A+|. 

(c) Spatial dependence of the local density of states for |A_| = 0.5|A + | in the unitary limit. 

Parameters are chosen as u++ = 2u = \/2u + _ and = 0.2£_. The ratio E F /A = 100 and the 

broadening factor T — 0.1A are used. 



function is denned by 

K(iv l )= f dr^K(r), K{r) = -(T t SL(t, 0)5+(0, 0)), (25) 
J o 

where v\ = 2irTl is a Matsubara frequency for boson. The spin operators are written by the 
field operators at the impurity site. 

S_(t,0)=^!(t,0)V T (t,0), S+(0,Q) =V}(0, 0)^(0,0). (26) 

For the multiband, the field operator is written as a summation of that for the \i = db bands. 42 ) 

W(0) = £ ^(0). (27) 

The two body Green's function is then written as 

K ( T ) = ~ E E (^^(^0)^^(^0)^(0,0)^^(0,0)} 

flfl'=± vv' = ± 

= ~ E E <rr^ 1 (r ) o)^i(o,o))(r r ^ T (r,o)^ T (o,o)) 
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+ E E (^(r,0)^ T (0,0))(T T ^ T (r,0)^ U (0,0)). (28) 

fj,fj,'=± w'=± 

Without the interband scattering, the expectation values are diagonal for the band index, for 
instance, (T T i/J ^(t, 0)Vvj.(0, 0)) = (T r -i//^(T, 0)Vyi(0, 0)) 6^. In contrast to this, off-diagonal 
elements remain in the presence of the interband scattering. Using the spectral representation, 
we can express T-f 1 at the impurity site as 

Tl tx y d£; i + cosh(£;/r) • (29) 

Here, g, g, /, and / are expressed as 

g(E) = a n {E, 0) + a 33 {E, 0) + a 13 (E, 0) + a 31 (E, 0), 

g(-E) = a 22 (-E, 0) + a 44 (-E, 0) + a 24 (-E, 0) + a 42 (-E, 0), 

f(E) = a 12 {E, 0) + a 34 (E, 0) + a u (E, 0) + a 32 (E, 0), 

f(-E) = a 21 (-E, 0) + a 43 {-E, 0) + a 23 (- J E, 0) + a 41 (-E, 0). (30) 

a mn (E,r) is defined by eq. (15). In a single band case, eq. (29) is expressed by only the first 
terms in the right hand side of eq. (30). In a pure (no impurity) multiband case, the second 
terms also remain and the Hebel-Slichter peak is suppressed due to the cancelation of the s±- 
wave order parameters. For the single impurity, we need the third and fourth terms (interband 
spin correlations) as well in the presence of the interband scattering in the same manner as 
the local density of states. 

We show the temperature dependence of T^ 1 for the identical multiband in Fig. 5. For 
u = 0, Tj -1 shows a small Hebel-Slichter peak just below T c due to the canceration of the 
coherence factor for the s±-wave state. T^ 1 is reduced with the increase of u, however, a peak 
appears at lower temperatures. This does not originate from the Hebel Slichter peak but does 
from the impurity-induced local boundstates, since the nuclear magnetic relaxation is possible 
via the local boundstates. The temperature at the peak position is related to the energy of the 
boundstates. 40 ^ For larger u, the boundstate energies decrease as in Fig. 3(a) and the peak 
position shifts towards the low temperature region as in Fig. 5. At the impurity site, impurity 
effects appear in the local density states strongly. It reflects in the T^ 1 considerably. Thus, 
the impurity site NMR is sensitive to the existence of the low energy boundstates and it is 
one of a suitable probes for unconventional superconductivity. 

3. Effective Single Band Model for Identical Multiband 

In this section, we focus on the identical multiband case and discuss why the low-energy 
states appear by the interband scattering on the basis of an effective single band model. This 
model enables us understand the essential role of the interband scattering for the s±-wave 
state. 
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Fig. 5. (Color online) Temperature dependence of for the identical multiband for various values 

of u. The scatterings are symmetric (u ++ = u = u_| = u). We use the temperature dependent 

order parameter obtained by solving the BCS gap equation. The broadening factor T = 0.1A(0) 
is used, where A(0) represents the order parameter at T = 0. 



3.1 Nonmagnetic scattering 

For the s±-wave, the Green's function and scattering matrix are defined by eqs. (3) and 
(10) in the 4x4 matrix form, respectively. We first diagonalize the scattering matrix. 



XT 



eff 



A^IJA 



( 



\ 



+ ^4 






u. 



where A is defined by 



A 
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V2 




















-H 




U+- 





















u ++ - 


U+- 



















-U++ + U. 
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-1 





\ 
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1 








1 


1 












V 








1 


- 1 J 







\ 



We next transform the Green's function by the matrix A 
Gf(ioj l ,k) = A- l G (iio l ,k)A 

( 

1 



+ + 







-A 





iwz + e fc 
A 





A 

i^l - e k 




-A 





(31) 



(32) 



(33) 



^1 - e k J 

We notice that G c ^ has the same form of a single band s-wave Green's function in a 4 x 4 
matrix form. Thus, the problem reduces to a 2 x 2 matrix form even if there is an interband 

scattering. In the reduced matrix form eq. (31), we notice that J7_| works as an Ising spin. 

Therefore, the interband nonmagnetic scattering plays the role of a magnetic scattering and 
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is pair breaking for the s±-wave superconductivity. 

For the symmetric scattering = U-\ ), one of the effective potential becomes zero, 

while the other is ±2C/_| as in eq. (31). This means that only one of the conduction electron 

forming a Cooper pair is scattered by the potential, while the other is not. Using the effective 
Green's function reduced to the 2x2 matrix form, we can obtain the same boundstate energy 
Eb defined in eq. (20). 

3.2 Classical magnetic scattering 

Let us consider here magnetic scattering of Ising type. The matrix U in eq. (10) has the 
following form: 



/ 



U 



\ 



T z 

J ++ 





T z 

J +- 








J ++ 





T z 

J +- 


J +- 





T z 

J ++ 








T z 

J +- 





T z 

J ++ 



(34) 



Here, J+ + and J+_ represent coupling constants of the intraband and interband scatterings, 
respectively. As in the nonmagnetic case, we obtain the following boundstate energies: 



E = ± 



1 + (f ++ 

where j ++ and j-\ are defined by 

J++ = *N+ J++, 



•i 2 
J+- 



(35) 



j+. = iry/N+Nlj^. (36) 

Comparing the boundstate energy for the nonmagnetic [eq. (20)] and Ising [eq. (35)] cases, 
we notice that they are equivalent under the following transformations: 

u++ < — u+-< — >j ++ . (37) 

This result implies that the roles of the magnetic and nonmagnetic scatterings are interchanged 
for the interband scattering in s-t-wave superconductors. This property has been reported 
by Golubov and Mazin who studied reduction of transition temperature of s-t-wave pairing 
superconductors. 43 ^ 

We elucidate this point by mapping the s-^-wave model to an effective single band one. 
As in the nonmagnetic case, U is transformed as 



u eS = A' l UA = 



'+- 



\ 



+ J. 






++ 





_ TZ i TZ 

J +- + J ++ 









TZ _| TZ 

J +- + J ++ 










\ 



(38) 



Since the 11 and 44 components of U eS are coupled via the order parameter terms in eq. (33), 
the interband magnetic scattering J^„ works as a nonmagnetic scattering in the effective 
single band model in the reduced 2x2 matrix form. Thus, the roles of the magnetic and 
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nonmagnetic scatterings are interchanged for the interband scattering in the s±-wave state. 



3.3 Relation between s±-wave and d x i_ y i ± id xy -wave superconducting states 

For the s±-wave state, there are two superconducting conduction bands with isotropic 
s-wave order parameters. The characteristic point is that the signs of the order parameters 
are opposite. It makes difference between the s-t-wave and the conventional s-wave states 
as we showed in the single impurity problem. This indicates that the s±-wave state has 
unconventional pairing nature. Since the s±-wave state is a fully gapped singlet pairing state, 
we focus on a d x 2_ y 2 ± id xy -w&ve state in this subsection and discuss the single impurity 
problem in a different point of view. 

The Hamiltonian eq. (1) is rewritten in the momentum space. 

Ti = Wkin + + Wimp, 

Ti-kin = ^ 6 fc C jfco- C Mfco-' 

H=± ka 

Wa = - Yl { C UA,-H + Cn-kic^n) , 

At=± fc 

n imp = Yj EE^' c ^Vfev. ( 39 ) 

where c^ ka and c^o are creation and annihilation operators for the conduction electron with 
momentum k and spin a for the [i = ± band. The Hamiltonian consists of three terms: Wkim 
Ha, and TL\ mp are for the kinetic energy, for the pairing interaction, and for the interaction 
between the conduction electron and the impurity, respectively. Since there is a rotational sym- 
metry around the single impurity, it is convenient to use the polar coordinate. We transform 
then the operator as 

C^a = yj^EH)^^' ( 4 °) 

where R represents the system size. I is the z component of the orbital angular momentum of 
the conduction electron, k is the wave number. 0^ is the angle from the wave vector measured 
from the fc^-axis. The Hamiltonian eq. (39) is then rewritten as 

jU=± k cr 

^^"EE^ ( c UoT c ifcO| + c^oic^ot) , (41) 
M=± k 

Wimp = ^ ^ ^ U W' C \kOa Vfe'Off- 

HH'=± fcfc' <r 

Here, is the Fermi wave number and the summation means Y2k = § Jo°° Wimp is 
composed of the operator for the I = angular momentum. Since the I = orbital is connected 
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to the impurity, we retain only the I = component in the Hamiltonian. 7i/^ represents that 
the total angular momentum of the Cooper pair is zero due to the s-wave nature of the 
s±-wave paring state. Next, we transform the operator by 



CfikOa 



V2 



(42) 



Here, c a k a and cp\. a are annihilation operators for fermion. The Hamiltonian eq. (41) is then 
written as 

7=a,/3 k cr 



(43) 



Ti 



imp 



* 

7,7'=a,/3 fcfc' °" 

Here, [7 77 / is given by 

U a j3 = Ufa 

To discuss physical meaning of eq. (43), we consider Hamiltonian 7i/\ for the d x 2_ y 2 + id xy 
wave pairing state. 



V 



\ / u. 
i 

~ 2 



/ 



V 



++ + U— + 2l7 4 
U++ + L7__ - 2l7 4 



(44) 



H A = - A fc (4 T c f _ fci + c_ fci c fcT ) . 



(45) 



Here, = Ae 2l< ^ is the momentum dependent order parameter for the d x 2_ y 2 + id xy -wave 
state. Substituting eq. (40) into eq. (45), we obtain 26 ) 



= 2 CTA (4far4fc,-<r + c 0fc,-aC 2 fc<x) 



(46) 



where irrelevant angular momentum components disconnected to the impurity are truncated 
here. The total angular momentum of the Cooper pair is expressed by +2, reflecting the chiral 
d x 2_ y 2 + id xy -w&ve character. We notice that eq. (46) has the same form of 7i/\ in eq. (43). 
Therefore, the a and (3 indices introduced in eq. (42) correspond to the I = and I = 2 
angular momentum in the d x 2_ y 2 + id xy -w&ve picture. Thus, the single impurity problem in 
the s±-wave state is equivalent to that in the d x 2_ y 2 ± id^-wave state. 



Let us consider the symmetric scattering case {U^ 



U = [7-1 = U) in the d x 2_ y 2 + 



id xy -w&ve picture. Since only Uqq = 2U is finite [U aa = 2U in eq. (44)], it can be mapped to 
a short-range scattering impurity problem in the d x 2_ y 2 ± id xy -w&ve state. It is known that 
the short-range scattering gives rise to a local boundstate of energy given in eq. (21). 39 > 40 ) It 
explains why the value of 2u appears in eq. (21). Therefore, it is natural to have the mid-gap 
state in the unitary limit for the symmetric scattering in s±-wave superconductors. 
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Fig. 6. (Color online) Local density of states at E — (mid-gap) in the unitary limit, (a) s±-wave 
case with the symmetric scattering, (b) For d x 2_ y 2 ± id^-wave. The impurity is located at the 
origin of the coordinate. The radius r is scaled by the coherence length £. The ratio E-p/A = 100 
and the broadening factor T — 0.08A are used. 



Although energy of the boundstate is same in the s-t-wave and d x 2_ y 2±id xy -w&ve states, a 
little difference between them appears in spatial dependence of the local boundstates around 
the impurity as shown in Fig. 6, since the a (/?) index introduced in eq. (42) is not the angular 
momentum / = (/ = 2). The real space Green's function for the d x 2_ y 2 + id xy -wave is given in 
the Appendix. The difference between the two cases can be seen only in a microscopic length 
scale (Fermi wave length). In a long length scale such as the superconducting coherence length, 
there is no significant difference between them. 

We mention here another different point between the s-t-wave and d x 2_ y 2 ± ici^-wave 
states. Since the d x 2_ y 2 ±id xy -w&ve state breaks the time reversal symmetry, electric current 
is induced by scatterings such as an impurity, surface, and domain wall. 44 ' 45 ) In contrast to 
this, the time reversal symmetry is not broken in the s±-wave state and such current is not 
induced. 

4. Quantum Spin and Kondo Effect 

It is known that magnetic impurities destroy the superconducting order parameter and 
suppress the superconducting transition temperature. Although these results are for conven- 
tional BCS superconductors, they hold also in the s±-wave superconductors when the magnetic 
scattering is intraband type. However, the interband type is open to further investigation. In 
this section, we examine effects of the interband magnetic scattering in s±-wave superconduc- 
tors in the identical multiband case using the Wilson's NRG method 14 ) which is reliable to 
study the Kondo effect also in superconductors. 15,16 ) 

In the same manner as the nonmagnetic scattering [Wi mp in eq. (41)], the Hamiltonian for 
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the magnetic impurity is expressed as 

Wimp = ^2 ^2^2 S ■ & <T<T ,J W ,c ikOa C IJ.'k'0<T', ( 47 ) 

fifJ.'=± kk' ua' 

where S represents the S = 1/2 spin operator for the impurity. J++ and J (J^ = J |_) 

are for the intraband (interband) magnetic scattering. 

We examine the magnetic impurity problem as in the d x 2_ y 2 + iJ^-wave case. 26 ) In Fig. 

7(a), for various values of 6 = /J++, we show the energy of the lowest-lying spin-singlet 

state with particle-hole degeneracy measured from that of the lowest-lying spin-doublet at low 

temperatures in respect of the relevant coupling J TC \ = (J++ + J^ ), where J++ = J and 

N + = N- = Nq are assumed. The meaning of J rc \ is described later. Although appearance 
of the boundstates can be understood qualitatively by the Ising spin case in §3.2, for the 
quantum spin, there are two energy scales characterizing the competition of superconducting 
pairing and Kondo-singlet formation. One is the superconducting energy gap A and the other 
is the Kondo temperature defined simply as Tk = Nq J++ exp(— 1/Nq J++) for 6 = in the 
unit of the half width of band. For 6 = 0, the ground state changes from the spin-doublet state 
to the spin-singlet as J++ increases. The Kondo singlet is stabilized for a large 7k /A only 

when J-| = 0. This resembles the case of a local 5 = 1/2 quantum spin in a conventional s- 

wave superconductor. Once J_| is finite, the spin-singlet energy merges into the spin-doublet 

one for a large J++. This implies that J_| destabilizes the Kondo singlet. Besides b ~ 0, 

the qualitative behavior is represented by the 6=1 case described by an s-wave scattering 
magnetic impurity coupled to the chiral d x 2_ y 2 ±id xy -wme superconductivity discussed below. 
On the other hand, for a small Tk/A, the doubly degenerate bound (spin-singlet) state appears 
in the superconducting energy gap, like a nonmagnetic impurity in §2. One can also find that 
Tk is estimated to be Nq J re \ exp(— NqJ tc \). 

Let us discuss the above result for the s±-wave state in terms of the d x 2_ y 2 + id xy -w&ve 
picture. 7ii mp can be mapped to the d x 2_ y 2 + id xy -w&ve model as in the nonmagnetic scattering 
given in eq. (43). There is the following relation in the coupling constants: 



( Joo ^ 

\ J20 = J()2 J 



J. 



++ 



[ 1 + b ) 

1 - b 

V / 



(48) 



Let us consider here that the scattering is only the interband type (J++ = J — ► 0). In this 

case (6 » 1), the coupling constants are Joo = — J22 in the d x 2_ y 2 + id xy -wave model [see eq. 
(48)]. This set of parameter means that one of the Joo and J22 is antiferromagnetic and the 
other is ferromagnetic. Figure 7(b) is a schematic of the groundstate in the d x 2_ y 2 + id xy -w&ve 
picture. We can see that there is no frustration and the spin doublet ground state is stabilized 
even for a strong coupling, since J22 < is ferromagnetic here. This picture still holds for 
J ++ = J / as long as J^ is finite. This means that the larger one of Joo or J22 is 
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Fig. 7. (Color online) (a) The coupling J^i dependence of the boundstate energy levels for S = 1/2 

local spin in the s±-wave state. The data are for b = J H /./++ (</++ = J ), the ratio of strength 

of the interband and intraband scatterings: 6 = (circle), 0.1 (square), 0.5 (diamond), 1.0 (up- 

triangle) and 2.0 (down-triangle). J rc i = (J++ + ) is the relevant coupling in the Kondo effect 

(see text) and J re i = 1.5iVoJ r ei is used for our NRG analysis, (b) Schematic of the groundstate in 
the d x 2_ y 2 + idxy-wave picture. Jqo is antiferromagnetic, while J22 is ferromagnetic. A represents 
the Cooper pairing interaction between the | (I = 0) and j (I = 2) conduction electrons. 



relevant. Therefore J Te \ discussed above is equivalent to Joo that stabilizes the Kondo singlet 
only with one orbital. 

5. Summary 

In this paper, we investigated single impurity effects in s±-wave superconductors. The 
main results of this paper are as follows: 

(1) Energy of the impurity-induced local boundstate is expressed explicitly as a function 
of strength of nonmagnetic interband and intraband scatterings [see eq. (16)]. The result for 
the single impurity problem is related to the energy level of the in-gap state that appears in 
the density of states for many impurities. 36 - 1 Spatial dependence of the local density of states 
provides information for a local probe such as scanning tunneling microscope. 

(2) We suggested impurity site NMR experiment as a powerful probe of the local bound- 
states induced by the nonmagnetic interband scattering. It would capture some features of 
the s±-wave superconductivity. 

(3) Roles of the magnetic and nonmagnetic interband scatterings are interchanged in the 
s±-wave superconductors. We elucidated this point by mapping the s±-wave multiband model 
to an effective s-wave single band one. 

(4) Appearance of the single-impurity-induced local boundstates in s±-wave supercon- 
ductors can be understood by a single impurity model in chiral fully gapped d x 2_ y 2 ± id xy - 
wave superconductors. The s±-wave pairing state has similar unconventional nature of the 
d x 2_ y 2 ± id X y-wave superconductivity. 
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(5) For a quantum magnetic impurity case, we found that the interband scattering desta- 
bilizes the Kondo singlet with two band electrons. Appearance of the boundstates can be 
understood by a classical spin qualitatively, while the boundstate energy depends on Ik/ A. 
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Appendix: Real space Green's function 

A.l s±-wave state 

In this appendix, we calculate the real space Green's function. We assume isotropic two 
dimensional conduction electron system. For the \x = ± band, the Green's function is given 
by 

G^r, 0) - -- ^ 2 + e 2 +A 2 e 

2vr J_ v 7_ oc uf + e 2 + A 2 
Here, r is the radius from the center position of the impurity, (ft k is the angle of the wave vector 
measured from the k x -axis. We divide eq. (A-l) into two parts. The first is proportional to 
\uj[ — A^pi, and the second is proportional to ep3. We perform the integral of these parts 
independently. The first is calculated as 

h = d(f> k de e V ^ * 

J-n J-oo + e 2 + A 2 

= J± f d fc P de f -J ^ eV^p+^-cos^ 

2^ M / J- 

= 27F I' 7 ''' 2 d fc f e i^ F +2m M in^r cos fe + fc2 F -2m M in Mi r cos </> fc \ 

7-7T/2 V / 

2vr 2 

[J (fc M +r) -fiFo^+r) + J (k^r) - \H (k^r)] , 



2^ M / 

where Jo(z) and Hq(z) are the 0th Bessel and Struve functions, respectively. They are defined 
by 

I f*/2 l /-tt/2 

J (z) = - d</> fc cos(zcos0 fc ), H (z) = - d<p k sm(zcos(j) k ). (A-3) 

In eq. (A-2), we used 



v/^+A 2 , ^± = ^ 2 F ±2m^. (A-4) 
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The integral of the second part is calculated as 



„ U, 2 + 6 2 + A 2 



C / 1 I \ i ^k^ F +2m M ercoscpj. 



de 7T— + 



I 2 = J* d<p k J°° de-„ % ^ e V^ +2m " ercos ^ 

= I f d(j) k [ ... 

2 J-TT J-OO V 6 - lJ V 6 + lJ V - 

2vr 2 

= i— [Jo(Vi-r) + i^o(^+ r ) ~ M k n- r ) + i^o(^- r )] ■ 
Using /i and ^2, we obtain the real space Green's function as 
G„(m,r,0) = *N ~ lUJl + AmPi [J (^+r) + iH (k li+ r) + J (***- r ) " ^o(V)l 

+ mN,^ [Jo(^+r) + i# ( - J (k^r) + itf (VO] ( A " 6 ) 
= ~ 7rN f^f~ l U J +( k »+ r ) + J-{k^r)\ [iu)i - A M pi) + [J+(/c M+ r) - J_(/c M _r)] i0^p 3 } , 
where J±(z) are defined by 

J±(z) = J (z)±iflo(z). (A-7) 

Equation (A-6) reduces to eq. (11) for r = 0, since Jo(0) = 1 and Hq(0) = 0. In the practical 
calculation, we perform the integrals in eqs. (A-2) and (A-5) numerically. It is convenient to 
introduce the following band dependent coherence length ^ and dimensionless radius f M : 

^ = 2a: = ^ a:' r = ^ (A - 8) 

We can rewrite k^±r in eq. (A-2) as 



A. 2 d x 2_ y 2 ± id xy - wave state 



For the d x 2_ y 2 + id xy -w&ve state, the order parameter depends on the wavevector (A& = 
Ae l2< ^). Green's function in a 2 x 2 matrix form is given by 
A ,. ^ 1 {uJl + e kh ~ Ae'^fcpi ik . r 

- r ^ r d ^+^-A^ ^ M , t (A . 10) 

wf + e z + A^ v 
The term proportional to A is different form the s±-wave case. This term is calculated as 

A p i20£ _ 

3 iy k^+2mercos4>j t 



/7T /*0O 
d(p k / d< 
-7T ■/ — C 



c^ 2 + e 2 + A 2 
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r/2 

2vr 2 A 



_ 27rA f I ^^2<j> k /' e i A /fe2+2min i rcos0 fc _|_ g-i-y/fep -2mifV cos </> fc \ 
7-7r/2 ^ ' 



2«j 



J 2 (fe+r) + i 



2H l (k + r) 



k+r 



J 2 (fc_r)-i 



H (k_r) 



Full Paper 

(A-ll) 

2H 1 (k^r) 
k^r 



where J 2 {z) and #1(2) are the second Bessel and the first Struve functions, respectively. In 
eq. (A-ll), we used 

^ujf + A 2 , k± = yjkl±2mini. (A- 12) 



Using Ii, I2, and I3, we obtain the real space Green's function for the d x 2_ y 2 ± id xy -wave. In 
the same manner as eq. (A-8), we introduce a coherence length for the d x 2_ y 2 ± id^-wave. We 
perform the integral in eq. (A-ll) numerically as in eqs. (A-2) and (A-5). 
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